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Abstract 

QCD predictions for moments of parton multiplicity distributions are discussed. The 
next-to-leading terms and conservation law give rise to the peculiar oscillating shape 
of some ratio of the moments. The similar shape has been found by moment analysis 
of hadron multiplicities. Experimental data, theoretical Monte Carlo models and phe- 
nomenological fits have been used for ee, hh, hA, AA reactions at high energies. 



In multiparticle production the primary observable is the muhiphcity distribution P„ which 
contains in the integrated form all the correlations of the system. First attempts to derive parton 
multiplicities in the double logarithmic approximation (DLA) of QCD provided extremely wide 
distributions compared to those found experimentally for hadrons. Later, the solution of the 
QCD equation for the generating function which takes into account higher order terms of the 
perturbative expansion and conservation law was found and more efficient way of studying 
the distribution shape based on moment analysis was proposed (see review [0). When applied 
to hadron multiplicities, it shows peculiar oscillating curves first derived for partons in QCD 0. 
Their stability for reactions initiated by different projectiles and targets is rather astonishing 

i- 

Let us describe briefly the theoretical flndings. The generating function is deflned by 
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SO that probabilities P„, factorial (Fg) and cumulant (Kg) moments are 
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where C,^ = n\/m\{n — m)! are the binomial coefficients. 

The equation for the generating function in gluodynamics |^ looks like an ordinary birth- 
death equation but with a typical singular QCD kernel 

— = 1^ dxKixh',[Giy + ln(l - + Inx) - (4) 

where 7q = 2Ncas/'n, is the coupling constant, A^c = 3 is the number of colours, K{x) = 
^ — (1 — x)[2 — x(l — x)], y = InQ'^/Ql, is the jet virtuality, (5o=const. The y-dependence 
enters G{z,y) through the energy dependence of P„. The variable z is omitted in eq.(^. 

For fixed coupling constant, this equation can be reduced to the algebraic one and, therefore, 
has an exact solution |^. For running coupling, one can use Taylor series expansion of the 
terms in the non-linear part of the integral. It is instructive to demonstrate it here because 
lowest order terms give rise to DLA, and higher orders provide modified leading logarithm and 
next-to-leading logarithm approximations (see [|l], 0). After the expansion of G at large y in 
non-linear terms is done and both sides of eq.(§) are divided by G{y), one differentiates both 



sides of @ and the non-linear integro-differential equation (^) is reduced to the differential 
equation 

(lnG')" = 7o[G-l-2/iiG" + /i2G"], (5) 

where hi = 11/24:, h2 ~ 0.216 are given by integrals of the kernel K{x) with corresponding 
weights. The various approximations are obtained by considering in the righthand side of (|^) 
only two (DLA), three (MLLA) or four (NLLA) terms. Higher order derivatives are omitted 
here. 

Using the relation and equating the terms of the same power of z on both sides of the 
equation, one gets from for the ratio of cumulant to factorial moments 

rr 7g[l-2M7 + /^2(gV + gy)] ... 



with 



7 ~ 70 - \hill + ^(4/^2 - hD^l (7) 



In DLA, one neglects all the terms of higher order in 70 and gets from (|^) 

fff = 1. (8) 

Let us note that such a monotonous behaviour of the ever positive Hg corresponds to the 
negative binomial distribution with its parameter k = 2. The experimental fits show much 
larger values of k and, consequently, more narrow distributions. 

For the realistic values of 7o=0.48, the ratio Hg as given by eqs. (^, exhibits the 
minimum at 

Qmin ~ ~ 5 (9) 

hilo 

and then increases. With higher order terms of Taylor series expansion taken into account the 
ratio Hg oscillates 0. 

These oscillations are due to oscillatory behaviour of cumulants which have the physical 
meaning of genuine correlations in the g-parton system. The location of the first minimum is 
determined by the condition 57 ^ 1 which is related to the well-known quantum-field theory 
problem of the break-up of the perturbative expansion due to the factorial increase of the num- 
ber of graphs for n-particle process. In traditional language, it corresponds to the expansion 
parameter a^ra appearing instead of as in that case. At the same time, this location is deter- 
mined by the singular nature of the kernel at x ^ (and x — > 1 in the symmetrized form). For 
theories with regular kernels like A0g Q the first minimum is located at larger values of q. 
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Figure 1: 



The KNO-scaling is valid in such an approach, in practice. The energy dependence of Hq 
should be extremely weak. It is only due to running property of 70 and almost cancels in the 
ratio Hq as seen from (P). 

Surely, there is no experimental information on parton distributions. The local parton- 
hadron duality is often used to relate them to hadron distributions by imposing proportionality 
assumption. Otherwise, the hadronization scheme should be developed. Anyway, one may not 
confront directly theoretical predictions about partons to experimental data. 

Nevertheless, QCD shows a way to the new method of analysis of experimental data. Cumu- 
lants are very sensitive to slight variations of the distributions because of subtractions implied 
in eq.(|^). Therefore, it is worthwhile to calculate factorial moments of experimental multiplic- 
ity distributions and then cumulants, using (^, to find out the behaviour of the ratio Hg for 
hadrons. Its values in e~^e~ collisions at Z° peak 10] and in hadronic reactions [|, ^ are 
shown in Figs.l and 2, correspondingly. 

One concludes that the oscillation pattern is observed in all the cases. The position of the 
first minimum is in the range g = 4 — 6. The value of Hq in it steadily increases (in modulus) 
from ee to pp, hA, AA. The model calculations according to the dual parton model (DPM) |rT| 
and to the quark-gluon string model (QGSM) |T^ were used for hA and AA in Fig. 2 but it is 
shown in @] that they fit experimental data rather well. On the contrary, no phenomenological 
fit (ranging from Poisson to modified negative binomial distribution) is able to reproduce such 
a pattern. They give rise either to monotonous positive ratios or to ones changing the sign at 
each value of the rank. 
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Figure 2: 

Thus, the moment analysis is a more powerful method than the direct fits of multiplicity 
distributions. The cumulants are very sensitive to any uncertainties and require high statistics 
experiments at high energies. Among different factors influencing the shapes of Hq-cmves we 
mention the cut-off of the high multiplicity tail of the distribution (which is mainly due to 
limited experimental statistics), difference between the values of moments for charged particles 
distributions and those for negatives only [Q, the similar difference (in model calculations) 
between clusters and their decay products, and, finally, the experimental selection criteria and 
error bars of P„. All of them disappear at asymptotically high energies and with statistics 
increased but they can be important for any given experiment. No careful analysis of all these 
factors has been done up to now. We rely on the qualitative similarity of Hq-cmves in all the 
above cases and suppose that it is related rather to the underlying dynamics of the processes 
than to varying from one experiment to another selection criteria and statistics. 

Let us mention at the very end another byproduct of the analysis of the generating functions. 
When the sum in eq. (|1|) is cut at some final multiplicity due to finite experimental statistics 
and conservation law, the truncated generating function becomes the polynomial of A^-th order 
in z, i.e. it possesses A^ complex conjugate zeros in z-plane. It happens that those zeros lie near 
the circle of the unit radius and at large A^ come very close to the real axis and to the point 
z = 1 (see 10]), where all the moments are calculated according to eq.(0). It explains why the 
moment analysis is so sensitive to tiny details of the distributions. Some statistical analogies 



arise also in connection with these zeros 13 
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